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Abstract. A new method is proposed to compute the bulk viscosity in strange quark 
matter at high densities. Using the method it is straightforward to prove that the bulk 
viscosity is positive definite, which is not so easy to accomplish in other approaches 
especially for multi-component fluids like strange quark matter with light up and down 
quarks and massive strange quarks. 
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There are a lot of works in literature concerning the bulk viscosity for quark and 
nuclear matter in compact stars. Some are based on the Urea and non-leptonic weak 
processes for nucleons and hyperons [IlElElElllEllBElElIIIlIiniinilTSlITlIISlIig 
[T71 [18] . Others are about Urea processes and d-s transitions for unpaired quark matter 
[la EQl EH ESI EHl El EH] and some color superconducting phases [26l EH EH]- For a 
brief review, see e.g. Ref. [29]. In the CFL phase the bulk viscosity comes from the 
lightest degrees of freedom, e.g. kaons, superfluid H modes etc., which are no longer 
quarks [301 EH [32]. In this note, we will only consider strange quark matter in normal 
state, for recent reviews of strange quark matter in stars, see e.g. [331 [M]. K is believed 
that in the regime of high density far above the saturation one of nuclear matter, the 
degrees of freedom become quarks and gluons. At weak coupling limit for quark matter 
one can make reliable predictions based on perturbative quantum chromodynamics. 
The bulk viscosity in compact star environment is partly due to damping mechanism 
for external perturbations which can lead to instabilities such as r-mode instability 
[351 [361 [371 [38l [39l [ID]. The stars can drastically spin down by gravitational radiation 
when the bulk viscosity is insufficient. 

Bulk viscosity arises from energy dissipation due to the change of fluid volume. 
The fluid will deviate from thermal and/or chemical equilibrium when it is compressed 
or rarefled. The physical pressure (i.e. the instantaneous one) will differ from its 
equilibrium value given by the equation of state at the same energy density. Energy 
has to be consumed to re-equihbrate the system. The bulk viscosity, denoted by is 
deflned through the rate of energy dissipation, 

(^diss) = - /V ■ ^?dt = -- r pvdt , (1) 

^ Jo T Jo 

where r = 2it/u) represents the pulsation period, V ■ v the divergence of fluid velocity, 
n the number density, P the physical pressure, and V = 1/n the speciflc volume. 

The bulk viscosity can be deflned in an alternative way in terms of complex 
variables. The deviation of the physical pressure, P, from the corresponding equilibrium 
one, Pcq, can be related to the bulk viscosity in the form, 

5P = P - Peq = -CcV • V , (2) 

where (c is complex. The physical bulk viscosity is then given by the real part 3?(Cc)- 

There are two characteristic time scales related to bulk viscosity, one corresponds 
to perturbation and the other to relaxation. The dominant contribution to dissipative 
energy loss comes from the re-equilibration whose time scale is comparable to that of 
external perturbations. For example, in strange quark matter which consists of up and 
down quarks which are massless, strange quarks which are massive and electrons, the 
time scale of the weak interaction in the system is comparable to the typical one for 
pulsations induced by external perturbations. To calculate the bulk viscosity, we can 
only consider weak processes. 

In literature about the bulk viscosity for quark or nuclear matter, the simple fluid 
picture is always implied. In the real world, however, the strange quark mass is not 
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negligible. For those external perturbations which do not distinguish particle species, 
the pressures felt by different ingredients are the same. This would lead to equal 
momentum variations but different velocity ones due to different masses. Although 
the time scale for strong interaction through which the quarks of different flavors reach 
thermal equilibrium is much shorter than that of weak interaction responsible for the 
bulk viscosity, there are no constraints about fluid velocities of different species either 
from hydrodynamic equations or from particle distribution functions, because the fluid 
velocity is the thermal average of particle microscopic velocities, which vary widely for 
particles with different masses even when they have the same temperature. Especially 
the divergences of the fluid velocities for particles with different masses are not necessary 
the same due to the difference in equations of state, since we have V ■ v = — from 
the fluid equation in the local co-moving frame, where the energy density e and the 
pressure P depend on masses. Therefore there is no strong reason that strange quark 
matter should be regarded as a simple fluid if the strange quark mass is large compared 
to temperature and chemical potential, i.e. the simple fluid picture for strange quark 
matter is just an approximation when the strange quark mass is very large. A natural 
question to ask is: whether conventional approaches in literature which are successfully 
applied to the simple fluid picture are still valid for multi-component fluids. 

Based on above arguments we extended conventional approaches to strange quark 
matter as a multi-component fluid. We obtained general constraints on fluid velocity 
divergences for different particle species enforced by baryon number conservation and 
electric charge neutrality. We also found a new oscillation pattern which gives a new 
solution to the bulk viscosity in the multi-component fluid satisfying these constraints. 
The resulting bulk viscosity is about one order of magnitude larger than those from 
conventional approaches. Our further investigations indicate that the positivity of the 
results in our previous paper [25] can not always be guaranteed for some cases although 
allowed by those constraints. In this note we propose an alternative way of computing 
the bulk viscosity of strange quark matter as a mult i- component fluid in terms of entropy 
production, where the bulk viscosity is definitely positive. The organization of the note 
is as follows. First we summarize the formalism to calculate the bulk viscosity for strange 
quark matter. Then we introduce our new approach to bulk viscosity by relating the 
energy dissipation to entropy production. 

In strange quark matter the weak processes via the exchange of bosons relevant 
to bulk viscosity are 

1. u + d ^ u + s , (d^s transition), 

2. u + e<-^d + i', (Ureal), 

3. M + e^s + z/, (Urea II). (3) 

Note that the reverse reactions in the second and third lines (Urea I and II) are not real 
ones but represent d/s ^ u + e + u. In thermal and chemical equilibrium all chemical 
potentials satisfy. 



fJ-d — fJ's , 
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jJ^u + jJ'C = fJ'S ■ (4) 

Considering there is no accumulation of neutrinos in the system, we have set fi,y/p = 
because neutrino mean free paths are much longer than typical radii of compact stars. 
The local charge neutrality also plays an important role, 

^ QiUi - rie = , (5) 

i=u,d,s 

where Qi are electric charges in unit of e. From Eqs. (jlj) and ([5]) the equilibrium values 
of chemical potentials and number densities can be obtained. 

Volume oscillations induced by perturbations lead to chemical non-equilibrium, 
since the chemical potential is related to the number density Ui for quark species i 
through the equation of state. The baryon number density is then ub = We 
can use the oscillation of the baryon number density, 6nB = (5nBoe"^*, to account for 
the volume oscillation. Hereafter we will keep the subscript i for quarks and j for 
quarks and electrons. Due to local charge neutrality, the density oscillations Snj are not 
independent from each other. The chemical potential differences between two sides of 
reactions are written as 

AyUl = yUs - /id = - Sfid , 

AyUs = /is - Ai« - Ate = Sfis - 5/i« - djJe ■ (6) 

The above equations can be put into more compact form, Afik = CkjSfij. Here Ckj 
[k = 1,2,3 denote different reactions) are coefficients in front of chemical potentials in 
(E]). The reaction rates are proportional to these differences of chemical potentials when 
A/Xfe <T(«:/i), 

Tk = Xk^fik (= XkCkjSfij) ■ (7) 
The coefficients can be obtained from Boltzmann equation with collision terms. 
For each ingredient in strange quark matter continuity equations read 

hj + Uj'V ■ Vj = Jj . (8) 
The source terms in the right-hand-side can be expressed in a linear combination of F^, 

Jj = —CkjTk = —CkjXkCkj'SfJ'j' ■ (9) 

Note that there is no source term for baryon number since it is conserved in all processes. 
So we have 

he + V • vb = . (10) 

The conservation of the baryon number and the electric charge in all three reactions 
lead to relations about sources, 

IX^- =0. 

i 
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Substituting above relations into continuity equations (jS]) and (fTUi) . we obtain general 
constraints for fluid velocity divergences, 

nijV ■ vb = ^ ^n^V ■ Vj , 

i 

neV-Ve =^(5i«iV-Vj. (11) 

i 

To solve Uj additional inputs are required. In traditional treatments ueXj = Jj is used 
implicitly, where Xj = rij/uB is the number density fraction of baryons for particle j. 
It can easily be verified from the general constraints ( ITTi) that the usage of ueXj = Jj 
implies V ■ Vj = V ■ for all particle species. This holds only if the mass differences 
among different particle species are negligible. However the strange quark mass nis could 
be much larger than light quark ones. So the fluid velocity of the strange component 
could differ from that of the light sector. In our previous paper [25] we considered a 
limit case where V ■ = and V ■ = V ■ v^, which provides additional conditions to 
flx Uj. The numerical results are about one order of magnitude larger than those from 
traditional solutions [IHl EQl EH ESI [231 Elj . 

Now we turn to the issue about the positivity of the bulk viscosity. We will 
provide an explicit and general demonstration about the positivity regardless of whether 
the system is a simple fluid or mult i- component fluid. The angle in which we look 
at the problem is the possible connection between the energy dissipation and the 
entropy production related to bulk viscosity arising from irreversibility of the periodic 
compression-expansion process. Note that in an irreversible process the entropy of the 
system always increases. 

To observe the entropy change we separate the fluid into small volume cells and 
focus on one of them. Since the fluid velocity divergences for particles with different 
masses are not identical, or the corresponding hydrodynamic motions are different, the 
cell is not a closed system and has exchange of particles with others. Due to baryon 
number conservation we can choose the cell in the way that the total baryon number 
in the cell does not change. It is an open system but with a global control number, 
so we can call it quasi-control number system. We assume that the quasi-equilibrium 
is reached in the cell, i.e. the system is in full thermodynamic equilibrium except 
the chemical one, since the chemical reactions from weak interaction proceed much 
slower than strong processes. The entropy of the cell increases during the chemical re- 
equilibration. Neglecting the exchange of particles and heat on the surface of the whole 
system with outside, the net change of particles and heat vanishes. Therefore after one 
cycle of compression-expansion process, the only thermodynamic quantity that changes 
is the entropy. We also neglect temperature change as in traditional treatments. The 
change of internal energy in one period reads 

/d f Td 

^t^t^Ts'j^N^^ = Nsj -^dt , (12) 

r ' T 

where the superscript / represents the local cell, A'^^ its baryon number, and sjy the 
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entropy per baryon number due to internal irreversible processes. Using the formula for 
entropy production in dissipative hydrodynamics we obtain, 

I^ = -a-^ = ^^- (13) 

Then the second type of definition for bulk viscosity becomes 

^ /(V ■ ^B?dt = [ —X,{Afi)ldt . (14) 
ns J J Ub 

T T 

Finally we get the positive bulk viscosity in the new approach, 
/.A,(A/i)idt 

^ - /,(V . .Bfdt ■ ^''^ 

It can be proved that the above formula can reproduce the bulk viscosity for the simple 
(single-component) fluids. 

In summary, from the connection between the dissipative energy loss and the 
entropy production a new approach to calculating bulk viscosity in strange quark matter 
has been developed. The positivity of the bulk viscosity can easily be observed in the 
approach for simple and mult i- component fluids. 
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